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Abstract. An optimal transport path may be viewed as a geodesic in the 
^ space of probability measures under a suitable family of metrics. This geodesic 

may exhibit a tree-shaped branching structure in many applications such as 
trees, blood vessels, draining and irrigation systems. Here, we extend the 
study of ramified optimal transportation between probability measures from 
Euclidean spaces to a geodesic metric space. We investigate the existence as 

I I well as the behavior of optimal transport paths under various properties of the 

metric such as completeness, doubling, or curvature upper boundedness. We 
also introduce the transport dimension of a probability measure on a complete 
geodesic metric space, and show that the transport dimension of a probability 
* measure is bounded above by the Minkowski dimension and below by the 

i-pH Hausdorff dimension of the measure. Moreover, we introduce a metric, called 

"the dimensional distance" , on the space of probability measures. This metric 

a gives a geometric meaning to the transport dimension: with respect to this 
metric, the transport dimension of a probability measure equals to the distance 
from it to any finite atomic probability measure. 

The optimal transportation problem aims at finding an optimal way to transport 
^ a given measure into another with the same mass. In contrast to the well-known 

^ Monge-Kantorovich problem (e.g. [T],[B],I7|,[Tl],[TS],[TH],I5D],[ll]),the ramified 

ly-j optimal transportation problem aims at modeling a branching transport network 

l/^ by an optimal transport path between two given probability measures. An essential 

I feature of such a transport path is to favor transportation in groups via a nonlin- 

ear (typically concave) cost function on mass. Transport networks with branching 
0^ structures are observable not only in nature as in trees, blood vessels, river channel 

networks, lightning, etc. but also in efficiently designed transport systems such 
^ as used in railway configurations and postage delivery networks. Several different 

approaches have been done on the ramified optimal transportation problem in Eu- 
clidean spaces, see for instance [16], [24], [19], [25], [26], g], [2], [27], [l3], [5], [28], 
^ and [29] . Related works on flat chains may be found in f23|, [T2|, f25| and [21] . 

This article aims at extending the study of ramified optimal transportation from 
Euclidean spaces to metric spaces. Such generalization is not only mathematically 
nature but also may be useful for considering specific examples of metric spaces 
later. By exploring various properties of the metric, we show that many results 
about ramified optimal transportation is not limited to Euclidean spaces, but can 
be extended to metric spaces with suitable properties on the metric. Some results 
that we prove in this article are summarized here: 
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When X is a geodesic metric space, we define a family of metrics da on the space 
A {X) of atomic probability measures on X for a (possibly negative) parameter 
a < 1. The space {A{X) ,da) is still a geodesic metric space when < a < 1. A 
geodesic, also called an optimal transport path, in this space is a weighted directed 
graph whose edges are geodesic segments. 

Moreover, when X is a geodesic metric space of curvature bounded above, we 
find in §2, a universal lower bound depending only on the parameter a for each 
comparison angle between edges of any optimal transport path. If in addition X 
is a doubling metric space, we show that the degree of any vertex of an optimal 
transport path in X is bounded above by a constant depending only on a and the 
doubling constant of X. On the other hand, we also provide a lower bound of the 
curvature of X by a quantity related to the degree of vertices. 

Furthermore, when X is a complete geodesic metric space, we consider optimal 
transportation between any two probability measures on X by considering the com- 
pletion of the metric space {A{X) ,c?q). A geodesic, if it exists, in the completed 
metric space is viewed as an a— optimal transport path between measures. The 
existence of an optimal transport path is closely related to the dimensional infor- 
mation of the measures. As a result, we consider the dimension of measures on X 
by introducing a new concept called the transport dimension of measures, which 
is analogous to the irrigational dimension of measures in Euclidean spaces studied 
by [13]. We show in §4.2.3 and 4.3.4 that the transport dimension of a measure is 
bounded below by its Hausdorff dimension and above by its Minkowski dimension. 
Furthermore, we show that the transport dimension has an interesting geometric 
meaning: under a metric (called the dimensional distance), the transport dimension 
of a probability measure equals to the distance from it to any atomic probability 
measure. 

In §5, when X is a compact geodesic doubling metric space with Assouad dimen- 
sion m and the parameter a > max |l — ^,0} , then we show that the space V {X) 
of probability measures on X with respect to da is a geodesic metric space. In 
other words, there exists an a-optimal transport path between any two probability 
measures on X. 

1. The da metrics on atomic probability measures on a metric space 

1.1. Transport paths between atomic measures. We first extend some basic 
concepts about transport paths between measures of equal mass as studied in 
with some necessary modifications, from Euclidean spaces to a metric space. 

Let {X, d) be a geodesic metric space. Recall that a (finite, positive) atomic 
measure on X is in the form of 



with distinct points Xi £ X and positive numbers rrii, where Sx denotes the Dirac 
mass located at the point x. The measure a is a probability measure if the mass 
^Jl^ = 1. Let A{X) be the space of all atomic probability measures on X. 

Definition 1.1.1. Given two atomic measures 



k 



(l.Ll) 




k 



(1.1.2) 
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on X of the same mass, a transport path from a. to h is a weighted directed acyclic 
graph G consisting of a vertex set V{G), a directed edge set E(G) and a weight 
function w : E{G) (0, +00) such that {xi, 2:2, 2^*;} U {j/i, j/2, J/n} C V{G) 
and for any vertex v € V{G), there is a balance equation 
(1.1.3) 

if V — Xi for some i — I, ...,k 



E 



w{e) — 'w{e) + < —rij, if v — yj for some j — 1, 



, n 



ei£E{G),e-=v eeE{G),e+=v 0, otherwise 

where each edge e E (G) is a geodesic segment in X from the starting endpoint 
e~ to the ending endpoint . 



Note that the balance equation (1.1.3) simply means the conservation of mass 
at each vertex. In terms of polyhedral chains, we simply have dG — b — a. 

Here, a directed graph G is called acyclic if it contains no directed cycles in 
the sense that for any vertex v ^ V (G) , there does not exist a list of vertices 
{vi,V2, ■ ■ ■ ,Vk} such that vi = Vk = v and [tij, tii+i] is a directed edge in E (G) for 
each i = 1,2, •••,/(;— 1. Reasons for introducing this constraint were given in |29l 
Remark 2.1.5]. 

For any two atomic measures a and b on X of equal mass, let Path{a, b) be the 
space of all transport paths from a to b. Now, we define the transport cost for each 
transport path as follows. 

Definition 1.1.2. For any real number a G (— oo, 1] and any transport path G G 
Path{a.,h), we define 

M„(G) w{e)''length{e). 

We now consider the following optimal transport problem: 

Problem 1. Given two atomic measures a and h of equal mass on a geodesic 
metric space X , find a minimizer of 

M„(G) 

among all transport paths G G Path (a, b). 

An Mq, minimizer in Path(aL, b) is called an a— optimal transport path from a 
to b. 

1.2. The da metrics. 

Definition 1.2.1. For any a G (— cx),l], we define 

da (a, b) = inf {M„ (G) : G G Path (a, b)} 

for any a, b ^A{X) . 

Remark 1.2.2. Let a and b be two atomic measures of equal mass A > 0, and let 
a = ^a and b = ^b be the normalization of a and b. Then, for any transport path 
Gg Path{a,h), we have G= {V {G) , E {G) , \w] is a transport path from a to 
b with Ma(G) = A"Ma(G). Thus, we also set da (a,b) = A^d^ (a,b) . 

It is easy to see that da is a metric on A{X) when < a < 1. But to show 
that da is still a metric when a < 0, we need some estimates on the lower bound 
of da (a, b) when a 7^ b. 
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We denote S{p,r) (and B{p,r), respectively) to be the sphere (and the closed 
ball, respectively) centered at p G X of radius r > 0. Note that for any transport 
path G, the restriction of G on any closed ball B(p,r()) gives a transport path 
^\b{p ro) between the restriction of measures. 

Lemma 1.2.3. Suppose a and b are two atomic measures on a geodesic metric 
space X of equal total mass, and G is a transport path from a to b. For each p X , 
if the intersection of G D S {p,r) as sets is nonempty for almost all r G [0,ro] for 
some To > 0, then 

pro 

(1.2.1) M„(G|b(p,,„)) > / Yl M^Td^, 

where for each r, the set 

:= {e eE{G) -.enSip^r)^ 0} 
is the family of all edges of G that intersects with the sphere S{p,r). 
Proof. For every edge e of G, let p* and p^ be the points on e such that 

d {p, p* ) = max {d(p,x) : X ^ e} and d{p,p^) = min {d (p, x) : a; G e} . 
Then, since e is a geodesic segment in X, 

length{e)>d{p*,p*)>\d{p,p*)-d{p,p^)\^ / XiAr)dr 



where Xi (*") is the characteristic function of the interval /g := [d(p,p*) ,d(p,p*)]. 
By assumption, E-r is nonempty for almost all r G [0, ro]. Also, observe that e G E^ 
if and only if Xi {'"') = 1- Therefore, 

Mo(G|b(p,.„)) = E [u'(e)]"length(e) 

> [w{e)f / Xi, {r)dr 

eeE(G\B{p,ro)) ° 
ro 

Y [w{e)]°'dr. 

eeEr 

□ 

The following corollary implies a positive lower bound on d^ (a, b) when a. ^ h. 
Corollary 1.2.4. Let the assumptions be as in Lemma \l.2.3\ and a <0. Then 

M„ (G|b(p,.„)) > A"ro, 

where A is an upper bound of the weight w (e) for every edge e in G\g(^p^^^y In 
particular, for any atomic measure 

k 

1=1 

on X with mass ||a|| := X^iLi ^ ^' ™^ have the following estimate 
(1.2.2) (a,||ap„) > ||a|r max {d{p,x,)} . 

\ ^ J l<i<k 
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Proof. When a < 0, we have for each r with nonempty, 

J2 [w(e)]">max[u;(e)]">A" 

where A is any upper bound of the weights of edges in G\gi^p^^y Therefore, by 

(imj, 

M„ (G|s(p,.„)) > A"ro. 

Now, let G e Path (^a, ||a|| S^^ . When a < 0, by the acychc property of transport 
paths, we have 

«;(e)< ||a|| 



for every edge e ^ E (G). Then, ( 1.2.2 1 follows by setting tq = niaxi<i<fe {d {p, Xi)} 
and A = ||a||. □ 



Lemma 1.2.3 also gives a lower bound estimate for positive a, which will be used 



m proposition 



3.0.151 



Corollary 1.2.5. Suppose < a < 1. For any a <eA{X) in the form of (1.1.1), 
p Cz X and ro > 0, we have 



(1.2.3) 



Proof Let A = E 



J2 rn, 

d{p,Xi)>ro 



< 



d{p,Xi)>ro 

for any < r < tq, we have 



Let G be any transport path from a to 5p. Then, 

eSE,. 

By lemma 1.2.3 since the function / (x) = x" is concave on [0, 1] when < a < 1, 
we have 



fro 

M„(G|s(p,.„)) > / M^Tdr 

•^0 eGE„ 



> 



E w{e) 



dr> / A"dr = A"ro. 



Therefore, we have (1.2.31. 



□ 



By means of corollary 1.2.4 the proof of fS^", Proposition 2.2.3] shows the fol- 
lowing proposition. 

Proposition 1.2.6. Suppose a < 1, and X is a geodesic metric space. Then 
da defined in definition 1.2.1 is a metric on the space A{X) of atomic probability 
measures on X. 

1.3. The da metric viewed as a metric induced by a quasimetric. When 
< a < 1, another approach of the metric da was introduced in [28] , which says 
that the metric da is the intrinsic metric on A{X) induced by a quasimetric Jq, . 
Let us briefly recall the definition of the quasimetric J„ here. 



function g : X X X — > [0, +oo) is a quasimetric on X if 9 satisfies all the conditions of a 
metric except that q satisfies a relaxed triangle inequality q (x, y) < C {q {x, z) + q (2, y)) for some 
C > 1, rather than the usual triangle inequality. 
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Let a and b be two fixed atomic probability measures in the form of (1.1.21 on 
a metric space X, a transport plan from a to b is an atomic probability measure 

m I 
1=1 j=l 

in the product space X x X such that 



(1.3.2) ^7,j=njand 



for each i and j. Let Plan (a, b) be the space of all transport plans from a to b. 

For any atomic probability measure 7 in X x X of the form (1.3.11 and any 
< a < 1, we define 



Ha (7) — XI 51 (7»j)" ^ i^i^Vj) ' 

where d is the given metric on X. 
Using Ha, we define 

Ja (a, b) := min {Ha (7) : 7 G Plan (a, b)} . 

For any given natural number G N , let An{X) be the space of all atomic 
probability measures 

rn 

1=1 

on X with m < N, and then A{X) ~ IJ^^Ar(X) is the space of all atomic 
probability measures on X. 

In [28l Proposition 4.2], we showed that Ja defines a quasimetric on An{X). 
Moreover, Ja is a complete quasimetric on An (X) if (X, d) is a complete metric 
space. The quasimetric Ja has a very nice property in the sense that this quasi- 
metric is able to induce an intrinsic metric on {X). 

Proposition 1.3.1. [551 Theorem 4.17, Corollary A:.lS\Suppose < a < 1, and X 
is a geodesic metric space. Then, the metric da defined in definition (1.2.1) is the 
intrinsic metric on An {X) induced by the quasimetric Ja. 

Moreover, in |2S1 remark 4.16] we have a simple formula for the cost. Suppose 
G £ Path {a, h) for some a, b gAn{X). If each edge of G is a geodesic curve 
between its endpoints in the geodesic metric space X, then there exists an associated 
piecewise metric Lipschitz curve g : [0, 1] An (X) such that 

Ma{G)= [ \g{t)\jyt. 
Jo 

where the quasimetric derivative 

\9(*)\.J.-=h] 

exists almost everywhere. 

Corollary 1.3.2. Suppose (X, d) is a complete geodesic metric space. Then, {An (X) , 
is a complete geodesic metric space for each < a < 1. 
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Since Ai {X) C A2 {X) C • • • C An (X) C • • • , and {An (X) , da) is a geodesic 
space for each N , we have the following existence result of optimal transport path: 

Proposition 1.3.3. [28l proposition 4.02] Suppose {X,d) is a complete geodesic 
metric space. Then, (A{X) ,da) is a geodesic metric space for each < a < I. 
Moreover, for any a, b €zA{X), every a— optimal transport path from a. to h is a 
geodesic from a. to h in the geodesic space {A{X) ,c?q). Vice versa, every geodesic 
from a. to h in (A (X) , da) is an a~ optimal transport path from a to b. 

2. Transportation in metric spaces with curvature bounded above 

In this section, we will show that when X is a geodesic metric space with cur- 
vature bounded above, then there exists a universal upper bound for the degree of 
every vertex of every optimal transport path on X. 

We now recall the definition of a space of bounded curvature [5]. For a real 
number k, the model space is the simply connected surface with constant 
curvature k. That is, if fc = 0, then is the Euclidean plane. If A: > 0, then 
is obtained from the sphere by multiplying the distance function by the constant 
If fc < 0, then is obtained from the hyperbolic space by multiplying 
the distance function by the constant The diameter of M| is denoted by 

Dk n/Vk for fc > and Dk '■— 00 for fc < 0. 

Let {X, d) be a geodesic metric space, and let A ABC be a geodesic triangle in 
X with geodesic segments as its sides. A comparison triangle AABC is a triangle 
in the model space M| such that d {A, B) = \A- B\^ , d {B, C) = |B - C]^ and 
d{A,C) = \A — C\i^, where denotes the distance function in the model space 
M^. Such a triangle is unique up to isometry. Also, the interior angle of AABC at 
B is called the comparison angle between A and C at B. 

A geodesic metric space {X, d) is a space of curvature bounded above by a real 
number k if for every geodesic triangle AABC in X and every point h in the 
geodesic segment 7^^, one has 

d{h,B) < \h- B\^ 

where h is the point on the side j^q of a comparison triangle AABC in such 
that \h-C\^ = d{h,C). 

Now, let X be a geodesic metric space with curvature bounded above by a real 
number k. Suppose a < 1 and C is an a— optimal transport path between two 
atomic probability measures a, b £A{X). We will show that the comparison angle 
of any two edges from a common vertex of G is bounded below by a universal 
constant depending only on a. Moveover, when X is in addition a doubling space, 
then the degree of any vertex w of G is bounded above by a constant depending 
only on a and the doubling constant of X. 

More precisely, let O be any vertex of G 
and ei be any two distinct directed edges 
with ef — O (or = O simultaneously) 
and weight > for j = 1, 2. Also, for 
i = 1,2, let Ai be the point on the edge 
Bi with d(0, Ai) — r for some r satisfying 
< r < \Dk and r < length (cj). 
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Now, we want to estimate the distance d (Ai, A2). To do it, we first denote 
(2.0.3) R 



(m" + 7712 )^ ~ ("^1 + ™2)^" 



if 1712 



(2.0.4) i?„ 



TO" 

and have the following estimates for R: 

Lemma 2.0.4. For each a < 1, the infimum of R is given by 

V2, z/0<a<i 
V4 - 4", i/ i < a < 1 or a < 0. 

For each < a < 1, the supremum of R is given by 

o V2, z/i<a<l 

" ■ \ V4^4^, i/O < a < i. 

^Zso, tii/ien a = 0, then R = \/3- When a = |, i/ien i? = ^/2. When a — 1, then 
R = 0. 

When a < 0, we will show R < 2 later in lemma [2.0.6| 
Proof. We first denote 

(2.0.5) = = 

TOi + m2 TOi + TO2 

as in [231 Example 2.1]. Note that fci + fc2 = 1 and 



By considering the function 

fa {x) = 



l.a i.a 



(a;" + (l-x)")^-l 



for X £ (0,1) and a < 1, we have R = fa (ki). Using Calculus, one may check 
that for each x £ (0, 1), 

(1) when a € (O, the function is strictly concave up and 

4 - 4" = (l) < fa {x) < lim fa {y) = 2; 

(2) when a £ (57 l)i the function fa is strictly concave down and 

4 - 4" = f ^"j > fa (x) > lim /„ (y) = 2; 

(3) when a G (—00, 0), the function is strictly concave down and 

fa{x)>fo. f^) =4-4"; 



2^ 

(4) /q, has constant values when a G {O, ^, l}. 

Using these facts, we get the estimates for R — ^ fa (k\) for each a. □ 

Now, we have the following key estimates for the distance d (v4i, A2): 

Lemma 2.0.5. Assume that d{0,Ai) ~ d{0,A2) = r and < r < \Dk. Then, 
we have the following estimates for a :— d{Ai, A2): 
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(1) lfk>0, then 



cos ( aV fc ) < 1 ^ sin^ | r 



(2) Ifk^O, then 

(3) // fc < 0, then 



Vk 
a > Rr. 



. a\/k ^ R . I I, 

I.e. sm > — sm rv fc 

2-2 



cosh 



> 1 



— sinh^ ( r\ 
2 



fc I i.e. sinh 



> — Sinn r\ 

2-2 



Proof. Let P be the point on the geodesic JA1A2 fro™ Ai to A2 with 

d{A,P) = KdiAi,B) 



for i — 1,2 and some Ai G (0, 1) to be chosen later in (2.0.7) with Ai + A2 = 1. For 
any t G [0, 1], let Q (t) be the point on the geodesic -fQp from O to P such that 

d{0,Q(t)) = tb and d{P,Q{t)) ^{l~t)b 

where h = d{0,P). For i = 1,2, let AAiPO be a comparison triangle of AAiPO 
in the model space M^. 




(a) A triangle AAiPO in X (b) A comparison triangle AAiPO in M| 

Figure 1: Comparison triangles 



Thus, 

\Ai - Aifl, |Ai - = r and |0 - P|^ = 5. 

Let Oi (t) be the point on the side jQp of AAiPO such that |0 — = 
d{0,Q (t)) = tb and let ai (<) — \Ai — Qi Since X has curvature bounded 

above by fc, we have cr^ (t) > d {Ai, Q{t)). Let 

H{t) : =- ^[m?ai(t) + m>2(<) + (mi+TO2)"ifo] 

(mi + TO2) 

= fc^ai (i) + fc>2 (t) + 

> 7 ^ [m'id {AuQ (t)) + m^d (A2, Q (t)) + (mi + m2)" d (O, g (t))] 

(mi + m2) 

m?d(0,Ai)+mgd(0,A2) 

> 7 ^ = (U) , 

(mi + m2j 
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since O is a vertex of an a— optimal transport path G. This iniphes that H' (0) > if 
H' (0) exists. Now, we may calculate the derivative H' (0) = fefcr'^ (0) + fc^cr2 (0) + 6 
as follows. 

When A: > 0, by applying the spherical law of cosines to triangles AA^PO and 
A^iQ (t) O, we have 

cos (^\ia^/lij — cos {r^/k^ cos {h^/k^ + sin (^^/k^ sin {b^/k^ cos 6*,; 

cos I^CTi (i) \/fc^ — cos |^r\/A:^ cos ^t6\/fc^ + sin ^r\/fc^ sin ^t6\/fc^ cos6'i, 

where 9i is the angle j^AiOP. Thus, 

sm cos {^ia^/k^ —sin ^6\/fc^ cos ^ct; (t) = — cos {r^/k^ sin ((1 -t) 6\/fc 

Taking derivative with respect to f at i = and using the fact cri (0) = r, we have 
cos {\ia\/k\ +sin (feVfc) sin {rVk\ a[ (0) A/fc = bVkcos (rVkj cos 



Therefore, for i = 1,2, 



6 cos ^r\/fc^ cos (^b^/ftj — 6 cos ^A,:a\/fc^ 
bVkj sin ^r\/fc^ 



Applying these expressions to H' (0) = fcf c'l (0) + fcj ctj (0) + 6 > 0, we have 
(2.0.6) 

(fc" + ) cos |^r\/fc^ cos |^6A/fc^+sin (^r\^^ sin > fc" cos (^XiaVk^ cos ^A2aA/A:^ 

By setting 

Ve'^' = (fc^* + k^) cos (r\/fc) + i sin (r\/fc) 
as a complex number, we have 

(fcj* + fc^) cos (^rVfc) cos (^foVfc) + sin ^rv^^ sin (^bVftj = V cos (Bi - bVk^ . 
On the other hand, as Ai + A2 = 1, we have 
k" cos (^XiaVltj + k2 cos ^A2a\/fc^ 
= fc" cos (^\ia\/~kj + fcj cos ^a\/fc^ cos {^Xia^fkj + ^in {^aVk^ sin ^AiaA/fc^ 
— + ^2 cos I^oa/A:^^ cos (^ia-\/k^ + fcj sin |^a\/fcj sin (^\la^/k 

= VFcos ( 62 — AioV^ 



where 



We 



»2 = ^fc^ + fco (^aVk^ ) + j (fc2 sin (a\/A;) ) = fc? + k^e' 



as a complex number for some 02 € [0,27r). Thus, inequality (2.0.61 becomes 
Vcos (Qi - bVPj > Wcos - Aia\/fc) . 
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Since < r < ^Dk, we have < a < 2r < n/Vk. Then it is easy to see that 
< 62 < aVk. Let 



(2.0.7) 



6? 

Ai = ^e(o,i) 



we have the inequality V > W. That is, 
(fcj* + k^f cos^ (^rVk^ + sin^ (^rVk^ > (k^ + cos {a^/k 
By simpHfying this incquahty, we get 



k'2 



" sin^ [a^/k^ . 



< 1 



■ 2 ( /t\ ,1 . aVk 
— sm rv fc and thus sm 



R . 

> — sm 
2-2 



The proof for the cases k — and fc < are similar when using the ordinary (or 
the hyperbolic) law of cosines in the model space M|. □ 



Using lemma 2.0.5 we have the following upper bounds for R defined as in ( 2.0.3 ) 
which is useful when a < 0. 



Lemma 2.0.6. Let R be defined as in (2.0.3). For any fc and a <\, we have 

R<2. 

Proof. By the triangle inequality, we have a < 2r < Dk. We now use the estimates 
in lemma [2.0.51 



When fc < 0, then 
sinh 

This yields i? < 2. 
When fc = 0, then 

so < 2. 

When fc > 0, then 



kr 1 > sinh 



— fca R 

> — smh 

2 - 2 



-kr 



sin [r^/k^ 



2r>a> Rr, 



. ayk R.I /- 
> sm — ^ — > — sm I rv fc 



as < ^ < r\/fc < f . Therefore, we still have R<2. 



□ 



The following proposition says that when a is negative, the weights on any 
two directed edges from a common vertex of an a— optimal transport path are 
comparable to each other. 

Proposition 2.0.7. If a < 0, then for each i ~ 1,2, 

fc, > ^ 



l + (l + 2")" 



where ki is defined as in (2.0.5). 

Proof. Without losing generality, we may assume that fc2 > fci. By proposition 



2.0.6 we have R<2. That is. 



(fcf + fc2°) 
La L,a 
"'I "-2 



1 



< 4. 
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Simplify it, we have 



^1 ^2 ' 



Since ki £ (0, ^] and a < 0, we have 



1 



< (fci)"" < 2". 



Siniphfy it again using fc2 = 1 ~ fci, we have 

ki > ^ 

l + (l + 2")~ 



□ 



We now may investigate the comparison angle 9 between Ai and A2 at O, given 
in figure \Ta\ 

Proposition 2.0.8. Let X be a geodesic metric space with curvature bounded above 
by a real number k. Let 9 be the comparison angle between Ai and A2 at O in the 
model space Af|. Then 



> arccos 1 — 



= arccos 



Thus, by {2.0. 4-), we have 



> 



1 1.2a h.2 

1. 2 



ifQ<a<\ 



arccos (2^"~^ if\<a<lora<Q 



Note that when = 0, this agrees with what we have found in j24l Example 2.1] 
for a "Y-shaped" path. Also, when a approaches —00, then 9a approaches tt, and 
when a approaches 1, then 9a approaches 0. 

Proof. When fc > 0, then by the spherical law of cosines. 



cos a 



cos 



sm I T 

2 



Vk 



< 



1 - ^ sin^ ( rVk 



(r^/k^ — cos^ (^r\/ltj 



'Vk 



2 ■ 



When fc < 0, then by the hyperbolic law of cosines 

— cosh (aV^fc) + cosh^ (^r^/—k^ 



cosf 



sinh (^r\/—k^ 



< 



— 1 — ^ sinh^ (r ' ' — "^'^ 



A:) + cosh (rV— fc) 



sinh {r\ 

When fc = 0, then by the law of cosines. 



-k) 



2 ' 



cos f = 



r'^ + r"^ — a? ^ 2r^ — R^r"^ 



2r2 



2r2 



= 1 - 



i?2 



□ 
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Now, we want to estimate the degree (i.e. the total number of edges) at each 
vertex of an optimal transport path. We first rewrite lemma 12.0.51 as follows. For 

2 

any real numbers a; < (^) and < y < 2, define 

4^ arcsin (| sin (y/x)) , if < x < (|)^ 

'ifix,y):^{ |, if.T = 

^sinh"^ (f sinh(y^)) if a; < 

Then, one may check that ^' is a continuous strictly decreasing function of the vari- 
able X and an increasing function of y. Moreover, for each fixed y, lYm.j;^_oo ^ (x, y) - 
1 and 

■^{x,y)> ^-ffj) ,2/ )= -arcsin ''^ 



^ ,2/ '"J TT V2/ 
By means of the function the lemma [2 . . 5| becomes 

Lemma 2.0.9. Assume that d{0,Ai) = d{0,A2) = r and < r < \Dk. Then, 
we have the following estimate for d {Ai, A2): 

2r>d (^1,^2) > 2r^ {r^k, R) > 2rCa 

where Ca ■— ^ arcsin (^) . 

Note that since lim.x^-ao {x, R) = 1, we have d{Ai,A2) is nearly 2r when k 
approaches —00. 
Let 

(2.0.8) ^{x,a) = l+ ^ ^J'"""^ 

for a; e M and a < 1. For each fixed a < 1, (k) ■— ^ {k, a) is a strictly increasing 
function of k with lower bound limfc^_oo {k, a) — 2, upper bound 



$(fc,a) < 1 



In 2 
and 



4'(0,a) = 1 



In 2 

As in [T71 10.13], a metric space X is called doubling if there is a constant > 1 
so that every subset of diameter r in X can be covered by at most Cd subsets of 
diameter at most |. Doubling spaces have the following covering property: there 
exists constants /3 > and Cp > 1 such that for every e e (0, every set of 
diameter r in X can be covered by at most Cpe^^ sets of diameter at most er. This 
function Cpe~^ is called a covering function of X. The infimum of all numbers 
/3 > such that a covering function can be found is called the Assouad dimension 
oi X. It is clear that subsets of doubling spaces are still doubling. For any subset 
K of X, let dim^i {K) denote the Assouad dimension of K. 

Theorem 2.0.10. Suppose X is a geodesic doubling metric space of curvature 
bounded above by a real number k. Let a < 1 and G be an a— optimal transport path 
between two atomic probability measures on X, and O is a vertex of G. Let deg (O) 
be the degree of the vertex O and r (O) be the maximum number r in (0, -^Dk] such 
that the truncated ball B(0,r)\{0} contains no vertices of G. Then, 
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(1) for any < r < r (O), we have 

deg(0)<2(C,)*(''''=^"), 

where Cd is the doubhng constant of X, and <i> is given in (2.0.1 

(2) Moreover, deg (O) < 2 (Cd)*^^'"'', which is a constant depends only on a 
and Cd- 

(3) If deg (O) > 2 (Cd)^, then the curvature upper bound 



lcg(0)/2> 

(4) In particular, if deg (O) = 2 (Cd)*^°'"\ then A: > 0. 

(5) If fc < 0, then 



■io)<\ 



dcg(0)/2 



Proof. For any < r < ^(O), let {A,} be the intersection points of the sphere 
S{0,r) in X with all edges of G that flows out of O. It is sufficient to show 

that the cardinality of {Ai} is bounded above by (Cd)*'-'^ By lemma 2.0.9 

{B [Ai,r^/ (r'^k,Ra))} are disjoint and contained in B (O, (l + * (r^fc, r) 
Since X is doubling, the cardinality of |i? (Ai,^' (r^k,Ra) f)} is bounded above 

by (Cd)*^''''''"). This proves (1). By setting r -> in (1), we have (2). Then (3) 
and (5) follow from (1), and (4) follows from (3). □ 

3. Optimal transport paths between arbitrary probability measures 

In this section, we consider optimal transport paths between two arbitrary prob- 
ability measures on a complete geodesic metric space {X, d) . Unlike what we did 
in Euclidean space [24], we will use a new approach by considering the completion 
of A{X) with respect to the metric da- Note that {A{X) ,da) is not necessarily 
complete for a < 1. So, we consider its completion as follows. 

Definition 3.0.11. For any a G (— oo, 1], letVa{X) be the completion of the metric 
space A{X) with respect to the metric d^. 

It is easy to check that (see O lemma 2.2.5]) if /? < a, then V^iX) C 
and for all /x, v in 'Pf){X) we have ci/3(/Lt, v) > da{fi, v). Note that when a = 1, the 
metric d\ is the usual Monge's distance on A^X') and V\ (X) is just the space V (X) 
of all probability measures on X. Therefore, each element in Va can be viewed as 
a probability measure on X when a < 1. 



By proposition 1.3.3 the concept of an a— optimal transport path on A{X^ 
coincides with the concept of geodesic in (^(X), da). This motivates us to introduce 
the following concept. 

Definition 3.0.12. For any two probability measures fi^ and /i^ on a complete 
geodesic metric space X and a <\, if there exists a geodesic in {Va {X) ,da) from 
/i^ to ii~ , then this geodesic is called an optimal transport path from /i+ to /i^. 

In other words, the existence of an a— optimal transport path is the same as 
the existence of a geodesic in ■Pc(X). Thus, an essential part in understanding the 
optimal transport problem becomes describing properties of elements of Va (X) , and 
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investigating the existence of geodesies in 7'a(X). Since completion of a geodesic 
metric space is still a geodesic space, by proposition 1 1 . 3 . 3] we have 

Proposition 3.0.13. Suppose X is a complete geodesic metric space. Then for 
any < a < 1, {Va {X) ,da) is a complete geodesic metric space. 

In other words, for any two probability measures fi^ , fjT € Va{X) with < a < 
1, there exists an optimal transport path (i.e. a geodesic) from /Lt"*" to In par- 
ticular, since atomic measures are contained in Va{X), there exists an a— optimal 
transport path from any probability measure fi € VaiX) to Sp for any p € X . 

A positive Borel measure fx on X is said to be concentrated on a Borel set A if 
li{X \ A) = 0. The following proposition says that if a is nonpositive, then any 
element of Va{X) must be bounded. 

Proposition 3.0.14. Suppose a < 0. If n (z VaiX), then fi is concentrated on the 
closed ball B {p, da (a*, Sp)) for any p Cz X . 

Proof. If /i e Va{X), then fi is represented by a Cauchy sequence {a„} e A{X) 
with respect to the metric da. For any p G X, by corollary |1.2.4[ each a„ is 
concentrated on the ball B {p,da {an,Sp)). Thus, /i is concentrated on the ball 
B {p,da{fJ.,Sp)) . □ 

When < a < 1, /X e VaiX) does not necessarily imply fi is concentrated on 
a bounded set. For instance, let X = M, and /i = ^S^n}, which is clearly 

unbounded. Then, da (/i,^{o}) = E^^i iJ2kLn ' ^ = E^^i (2^^)" < ^^'^ 
/i £ VaiX). Nevertheless, the following proposition says that the mass of every 
measure in VaiX) outside a ball decays to as the radius of the ball increases. 

Proposition 3.0.15. Suppose < a < 1 and fi — Xfi for some fi G Va (X) and 

A > 0. Then, for any point p £ X and r > 0, we have 

[MX\B(p,r))]" < M^iMp). 

In particular, if r > da{jj., XSp)^^" , we have 

(3.0.9) jliX\B{p,r)) <daijl,XSp). 

Proof. By corollary |1. 2. 5| 

[^,(X\B{p,r))]''<^^^^. 

Now, for any A > 0, 

[A(X\B(p,r))]" = [X^iiX\B{p,r))]" < x^'^dtM ^ MMM. 

□ 

4. Transport Dimension of measures on a metric space 

Now, a natural question is to describe properties of measures that lie in the 
space Va{X). The answer to this question crucially related to the dimensional 
information of the measure /z. So, we study the dimension of measures in this 
section. To do it, we first study properties of measures belong to a special subset 
of Va {X), and then using it to define the transport dimension of measures. The 
work in this section generalize results of from Euclidean spaces to a complete 
geodesic metric space, while the study of [53] is motivated by the work of [TB] . 
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4.1. (ia admissible Cauchy sequence. 

Definition 4.1.1. Suppose X is a complete geodesic metric space. Let {afc}^j^6e 
a sequence of atomic measures on X of equal total mass in the form of 

i=l 

for each k, and a < 1. We say that this sequence is a da — admissible Cauchy 
sequence if for any e > 0, there exists an N such that for all n > k > N there 
exists a partition of 

a -Va'^'^ 

1=1 

with respect to a^ as sums of disjoint atomic measures and a path (see figure^ 

Gl,ePath{m\'^6,^^.^,ai'}) 
for each i = 1, 2, • • • , N^. such that 

i=l 

Also, we denote Gf^ = i' which is a path from a^ to a„ with Mq,(G^) < e. 

Each da— admissible Cauchy sequence corresponds to an elem,ent in Va{X). 
Let 

VaiX) C VaiX) 

be the set of all probability measures fj, which corresponds to a da admissible Cauchy 
sequence of probability measures. For simplicity, we may write VaiX) as Va. 



Figure 2: An example of a transport path between and 




It is easy to see that if for each k, there is a partition of a^+i = J2i^=i ^i+i i 
th respect to a^. as sums of disjoint atomic measures and a path G^_|_2 ^ G 

Path (^mf'^S^(k-) , a^*^j^ ^ for each i = 1, 2, • • • , Nk such that 

'Nk \ 




< +00, 
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then {a„} is a da-admissible Cauchy sequence. 

Also, note that if /i, G I?q,(X), one automatically has da (a*, i^) < oo. 

4.2. Relation with Hausdorff dimension of measures. Let W denote s di- 
mensional HausdorfF measure on X for each s > 0. By means of corollary 1.2.4 and 
(3.0.91, the proof of 29, theorem 3.2.1] is also valid for the following theorem: 

Theorem 4.2.1. Suppose X is a complete geodesic metric space. If ^ G Pq, (X) for 

some a G (— oo, 1), then fj, is concentrated on a subset A of X with TL^-° [A] = 0. 

Definition 4.2.2. For any probability measure ^ on a complete geodesic metric 
space X , the Hausdorff dimension of fi is defined to be 

di-mn (m) = inf {dimn (A) : fi iX\A) = 0} , 

where dim^f (A) is the Hausdorff dimension of a set A. 

Thus, by Theorem |4.2.1| we have 

Corollary 4.2.3. Suppose X is a complete geodesic metric space. For any a < I 
and any fi G T)a{X), we have 

dim/f(/i) < — — . 

1 — a 

4.3. Minkowski dimension of measures. A nested collection 

(4.3.1) ^ =^{Q1 : i = 1,2,-- - ,iV„ and n = 1,2,---} 

of cubes in X is a collection of Borel subsets of X with the following properties: 

(1) for each Q", its diameter 

(4.3.2) Cict" < diam {Q^ < Csct" 

for some constants C2 > Ci > and some a G (0, 1). 

(2) for any k, I, i,j with I > k, either Q,^ n Q^- = or C Q^.; 

(3) for each Q"^^ there exists exactly one Q" (parent of Q"^^) such that 

-in+l 



(4) for each there exists at least one QJ+^ (child of Q^) such that g"+^ 



C 



Q" 

Each Q" is called a cube of generation n in T. If two different cubes Qf and Q" 
of generation n have the same parent, then they are called brothers to each other. 

In next section, we will see that for each bounded subset of a complete geodesic 
doubling metric space, there always exists a nested collection of cubes which covers 
the set. 

Definition 4.3.1. For any nested collection J- , we define its Minkowski dimension 
(4.3.3) dimM (T) := lim J^^^ 

n^oo log (^) 

provided the limit exists, where Nn is the total number of cubes of generation n. 

Definition 4.3.2. A Radon measure fj, on X is said to be concentrated on a nested 
collection T in X if for each n, 
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Definition 4.3.3. For any Radon measure n, we define the Minkowski dimension 
of the measure fj, to be 

dim TV/ in) := inf {diniAf (T)} 
where the infimum is over all nested collection T that ji is concentrated on. 

The proof of [29, theorem 3.3.5] is still valid for the following theorem: 

Theorem 4.3.4. Suppose fi is a probability measure on a complete geodesic metric 
space {X,d). //dimjv/(M) < jz^ for some < a < 1, then fi € {^)- 

4.4. Evenly concentrated measures. Now, we aim at achieving a similar result 
as in theorem |4.3.4| for the case a < 0. To do it, we introduce the following 
definition: 

Definition 4.4.1. Let (X, d) be a complete geodesic metric space. A Radon mea- 
sure fi on X is evenly concentrated on a nested collection T in X if for each cube 
Q" of generation n in T, either Q" has no brothers or ii{Q^^) > for some 
constant A > 0. 

Here, Q" has no brothers means that the parent of Qf has only one child, namely 
Qf itself. 

Some examples of evenly concentrated measures have been given in |29| . In 
particular, if /x is an Ahlfors regular measure concentrated on a nested collection 
in X, then /i is evenly concentrated on J^. 

Definition 4.4.2. For any Radon measure /i, we define 

dimu := inf {diniM (^)} 

where the infimum is over all nested collection T that \i is evenly concentrated on. 

Obviously, 

diniM (m) < dim;7 (/i) . 
The proof of [29l theorem 3.4.6] is also valid for the following theorem: 

Theorem 4.4.3. Let {X,d) be a complete geodesic metric space. Suppose ^ is a 
probability measure with dim[/(/x) < for some a < 1, then fj, € Da i^)- 

4.5. Transport dimension of measures. We now introduce the following con- 
cept: 

Definition 4.5.1. Suppose X is a complete geodesic metric space. For any proba- 
bility measure fi on X , we define the transport dimension of fi to be 

dimy (m) := inf i : M ^ T^a(X) 

a<l [_ 1 — a 

Note that if > dim^ (/i), then /i G and thus da {fj,,So) < +oo for 

any fixed point O E X. If in addition a > 0, then there exists an a— optimal 
transport path from ^ to So- 

By theorems 4.3.4 4.4.3 and 4.2.1 we have (see [55| theorem 3.5.2]) 

Theorem 4.5.2. Suppose X is a complete geodesic metric space. Let fi be any 
probability measure on X, then 

dim/f (p) < diniT (/i) < max{dimA/(^), 1}. 
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Moreover, we also have 

dm\H (/i) < dini7'(/i) < dini[/(A*)- 
In Example 3.5.3], we showed that for the Cantor measure /z, we have 

In 2 

dim_ff (^) = dimT(Ai) = AYmu{^) = . 

in 3 

4.6. The Dimensional Distance between probability measures. In this sub- 
setion, we will give a geometric meaning to the transport dimension of measures. 
Let {X,d) be a complete geodesic metric space. For any a < 1, let 

5„ {X) ^{K{^l~v)■.K>Q,^l,ve v^{x)} 

be a collection of signed measures. Clearly Sa^ {X) C Sa^ {X) if ax <a.2- 

Definition 4.6.1. Let (X, d) he a complete geodesic metric space. For any two 
probability measures fi, v on X , we define 

Din,iy) := inf{-^ : fi - u e Sa (X)}. 

a<l 1 — a 

The proof of [29j proposition 4.05] still valid for the following proposition 

Proposition 4.6.2. Let (X, d) be a complete geodesic metric space. Then, D is a 
pseudometri^on the space of probability measures on X. 

In general, D is not necessarily a metric. Indeed, for any two atomic probability 
measures a, b, we have a — b S 5^ {X) for any a < 1. Thus, D (a, b) = while a 
and b are not necessarily the same measure. Nevertheless, we may easily extend 
the pseudometric Z? to a metric on equivalent classes of measures. To this end we 
define a notion of the equivalent class on measures. 

Definition 4.6.3. For any two probability measures fi and v on X , we say 

^ ^ V if D{ii, v) — 0. 
The equivalent class of fi is denoted by [fi\. 

For instance, all atomic probability measures are equivalent to each other. 
Definition 4.6.4. For any equivalent class [fi] and [i>], define 



From this definition and proposition |4 . 6 . 2| clearly, we have the following theorem. 

Theorem 4.6.5. Let {X, d) be a complete geodesic metric space. Then, T) is a 
metric on V {X) / ^ . 

Definition 4.6.6. The metric D is called the dimensional distance on the space 
V (X) / ~ of equivalent classes of probability measures on X. 

We now give a geometric meaning to transport dimension of measures. 

Theorem 4.6.7. Let {X, d) be a complete geodesic metric space. For any positive 
probability measure ^, we have 

dimT(Ai)=D([^],[a]) = i?(Ai,a) 

where a is any atomic probability measure. 



pseudometric D means that it is nonnegative, symmetric, satisfies the triangle inequality, 
and D (fi, p.) = 0. But D (p, v) = does not imply n = i/. 
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Proof. Since a gPq i^) for any a < 1, we have 

D{^,a) = inf{— : ^-ae 
a<l I — a 

= inf{-^ :^eI?„(X)}=dimT(At). 
a<l 1 — a 

□ 

This theorem says that the transport dimension of a probability measure fi is the 
distance from fi to any atomic measure with respect to the dimensional distance. 
In other words, the dimension information of a measure tells us quantitatively how 
far the measure is from being an atomic measure. 

5. Measures on a complete doubling metric space 

In this section, we will show that there exists an a— optimal transport path 
between any two probability measures on a compact doubling geodesic metric space 
X whenever max |l— ;^,0}<a<l, where m is the Assouad dimension of X. 

Recall that a space of homogeneous type ([H]) is a quasimetric space X equipped 
with a doubling measure v, which is a Radon measure on X satisfying 

v{B{x,2r)) < Cv{B {x,r)) 

for any ball B {x, r) in X and for some constant C > 0. When {X, d) is a metric 
space equipped with a doubling measure ly, then the triple {X, d, v) is called a 
metric measures space. Recently, many works (see [5], [T7], etc) have been done on 
studying analysis on metric measure spaces, in particular, when the measure /i is 
doubling and satisfying the Poincar inequality. 

In [9] and [10], Christ introduced a decomposition of a space of homogeneous 
type as cubes and proved the following proposition: 

Proposition 5.0.8. Suppose (X, v) is a space of homogeneous type. For any A; G Z, 
there exists a set, at most countable Ik and a family of subsets QgCX with 9 G Ik, 
such that 

(1) ly {X \ UgQ^) = 0, Vfc e Z; 

(2) for any k, I, 9, 77 with I < k, either n Qj, = or Q^' C Qj^; 

(3) for each Q^~^^ there exists exactly one Qg (parent of Q^^^) such that 

(4) for each Qg there exists at least one Q^^^ (child of Qg) such that Q^^^ C 

Qg'i 

These open subsets of the kind Qg are called dyadic cubes of generation k due 
to the analogous between them and the standard Euclidean dyadic cubes. A useful 
property regarding such dyadic cubes is: there exists a point Xg £ X for each cube 
Qg such that 

S(4,Coa'=) CQ^ CB(4,Cia^) 
for some constants Co, Ci and a G (0, 1). Moreover, for any Xg and x'^, d (xg, x^) > 

Since z/ is a doubling measure, from (1), we see that X = UgQg, where Qg 
denotes the closure of Qg. Then, it is easy to see that there exists a family of Borel 
subsets {Bg] with 9 G Ik such that Qg Q Bg C Q'^ with X — UgBg for each fc, 
and {Bg} still satisfy conditions (2,3,4) above. 
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A very useful fact is pointed out in |17[ theorem 13.3]: every complete doubling 
metric space {X, d) has a nontrivial doubhng measure on it. Thus, one may also 
construct a family of disjoint Borel subsets {Bg} for {X,d) as above. Now, for 
any bounded subset X of X, we set to be the collection of all dyadic cubes Bg 
that has a nonempty intersection with K. It is easy to check that J-'k is a nested 



collection of cubes as defined in (4.3.1 1. Moreover, for any /? > dim^i (K), from the 



definition of Assouad dimension, we see that the cardinality Nn of all dyadic cubes 
of generation n that intersect 
some constant Cp > 0. Thus, 



of generation n that intersect with the set K is bounded above by Cp (cr") ^ for 



dimM^<lim ^°g^^(f)"' =/3. 

log^ 

This shows that dim7\/ J-k < dim^i (K) . 

Proposition 5.0.9. Suppose (X, d) is a complete geodesic doubling metric space. 
If fx is a probability measure concentrated on a bounded subset K of X , then 

diniji/ (/i) < dim^ {K) . 

If in addition, fi is Ahlfors regular, then 

dim[/ (/i) < dimyi (K) . 

Proof. Since /i is concentrated on K, we have /i is concentrated on the associated 
nested collection J^k- Thus, 

diniAf (/i) < dimM J^k < dim^ (K) . 

When fi is Ahlfors regular, /i is evenly concentrated on J-'k, thus dim;/ (/i) < 
dimAiK). □ 

In particular, we have 

Theorem 5.0.10. Suppose X is a complete geodesic doubling metric space with 
Assouad dimension m, and /i is any probability measure on X with a compact 
support. Let ^ — ^ < a < 1. Then, 

(1) /i G Va (X) if a > 0. In particular, if in addition X is compact, then 

{X)=Va{X)^ViX). 

(2) fi e Va (X) if II is Ahlfors regular. 

Proof. Let K be the support of /i. Then, dim^,/ (fJ-) < dim^ (K) < dim^ (X) — m. 
By theorem |4.5.2[ dimy (/i) < max {1, dim^f (/i)} < max{l,m}. Therefore, for 
any maxjl — ;^,0} < a < 1, we have > max{l,m} > dimT{fi), and thus 
/Lt € Va {X). When fi is Ahlfors regular on X, we have dim^ (/i) < dim^/ (/i) < 
dim^ (K) < m. Thus, if > m, then ^ e (X). □ 

Thus, by proposition |3. . 13) we have 

Corollary 5.0.11. Suppose X is a compact geodesic doubling metric space with 
Assouad dimension m. Then, the space {V{X),da) of probability measures on 
X is a complete geodesic metric space whenever maxjl— ;;^,0} < a < 1. In 
other words, there exists an a— optimal transport path between any two probability 
measures on X. 
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